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I re-examine mass estimates of the Coma cluster from pre- ASCA X-ray spectral 
observations. A large range of model dark matter distributions are examined, under 
the assumptions of hydrostatic equilibrium and spherical symmetry, to determine 
the widest possible allowed range on the total mass of the cluster. Within a 
radius of 1 Mpc, the total cluster mass is tightly constrained to be (6.2 ± 0.9) X 
10 14 Mq and the ratio of luminous baryonic matter to total matter lies between 
13% and 17%. Within a radius of 3 Mpc the total mass is (1.3 ± 0.5) X 10 15 M Q 
and the luminous matter fraction is 20%-40%. I find that tjm "universal" dark 
matter density profile proposed by Navarro, Frenk, and White^J based on N-body 
simulations of a standard cold-dark-matter dominated Universe, predicts a steep 
temperature gradient within the core of the cluster that is a poor fit to the Coma 
data and can be rejected at greater than 99% confidence. 



1 Introduction 

Determination of the total masses of galaxy clusters is of great interest because 
of what these measurements can tell us about the amount, nature, and distri- 
bution of dark matter. The ratio of luminous to dark matter, which provides 
a lower limit to the baryon fraction, appears to be considerably higher in rich 
clusters of galaxies than that allowed by Big-Bang nucleosynthesis in a flat 
(O = 1) A = Universe. The Coma Cluster, because it is nearby and well 
observed, stands as an important laboratory for studies of this kind. 

I was motivated by several factors to revisit measuring the mass of the 
Coma cluster from X-ray observations. One was to address the reliability and 
accuractof X-ray-derived cluster mass estimates, which have been questioned 
recentlyEl Another was to see if Coma's radial gas temperature profile, which 
is essential for determining the gravitating mass, showed evidence for a re- 
markably xapid drop at the cluster outskirts as seen in the cluster Abell 2163 
by j45CAlI The availability of previously unpublished scanning data from the 
Ginga satellite that provides strong constraints on the temperature distribution 
over spatial scales up to 1° was a further motivation. 

In this article I utiliae all aj-ailable pie-ASCA temperature data on Coma 
from Tenmatl EXOSATp GingaB and coded aperture mask data from Spacelab- 
IJI3 (sec Fig. |l|). The ROSAT surface brightness profile a was used to define 
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Figure 1: Projected temperatures of the hot gas in the Coma Cluster (top panel) from various 
sources with one-sigma error bars (Ginga scan data on the left part of the plot; Spacelab-2 
data on the right; and pointed data from Tenma, Ginga, and several EXOSA T fields in the 
middle). The temperature values from the overall best- fit dark matter profile are shown as 
the starred symbols. The bottom panel shows the residuals between the data and the model. 



the gas deputy profile and the galaxy distribution was taken from Millington 
and PeachE3 I determine the cluster mass by constraining the parameters of a 
dark matter halo function using methods described previouslyu Briefly, given 
the RO SAT gas density profile and a particular set of parameter values for the 
dark matter halo, one directly solves for the temperature profile, which is then 
compared to the observed temperature values using a figure-of-merit function 
(I employ x 2 here) to assess the quality of the fit. A value of H = 50 km s _1 
Mpc -1 , which yields a distance to Coma of 140 Mpc, is used throughout. 



2 Parameterized Dark Matter Halos 

For the dark matter profile I use the three-parameter functional form given by 
Pdm = /°o[l + {R/Rdm) 2 ]~ an d the allowed ranges of all three parameters 
are determined. The fits to the data are reasonably good (minimum x 2 = 21.4 
for 15 degrees of freedom; this model is plotted in Fig. ^) and are summarized 
in Fig. 0. The left panel shows that the allowed values of the scale length 
i?DM and index a are strongly correlated. Formally I can reject a dark matter 
model that is distributed like the X-ray gas with high significance (>99% 
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Figure 2: (Left) Allowed ranges (at the 90% and 99% confidence levels) of the scale length 
of the dark matter distribution (-Rdm) as a function of the power-law index of the matter 
distribution, a. (Right) Allowed ranges (also 90% and 99% confidence levels) of the masses 
of the luminous and dark matter in the Coma cluster within three fiducial radii as indicated 
plotted as a function of a. 



confidence), but a model with the dark matter distributed like the galaxies is 
allowed, although it is not highly favored. Note that the current data do not 
constrain the maximum value of a and, in fact, \ 2 continues to drop, albeit 
very slowly, for even larger values of a than shown. 

The total gravitating mass of Coma is well constrained by these fits as 
shown in the right panel of Fig. ^. Here the masses, integrated within three 
fiducial radii (1 Mpc, 3 Mpc, and 5 Mpc), were determined based on the 
allowed ranges of -Rdmj ol, and po- (Note how the dark matter masses reach 
asymptotic values at large values of a so, even though a is not constrained 
at the upper end, my mass estimates are robust.) The total cluster mass 
(luminous plus dark) as well as the baryon fraction within 1 Mpc and 3 Mpc 
are quoted in the abstract for the 99% confidence level. The mass within 5 
Mpc is (1.9 ± 0.9) x 10 15 M Q and the luminous matter fraction is 20%-55%. 

The temperature profiles we obtain are nearly all convectively stable, i.e., 
| — | < |, within the maximum observed extent of the cluster ~3.3 Mpc. 
The few models that violate this constraint lie between the 90% and 99% 
confidence contours at low values of Rdm and a and thus are statistically less 
likely to be acceptable solutions anyway. In summary, I find no evidence for 
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steep temperature gradients in Coma that might indicate nonhydrostatic or 
other exotic conditions. 



3 "Universal" Dark Matter Halos 

Based on their TV-body simulations xif a standard cold-dark-matter dominated 
Universe, Navarro, Frenk, & White til (NFW) find that the radial dark matter 
density profiles of systems ranging from dwarf galaxies to rich clusters of galax- 
ies can be well described by the simple function p = S c p CI it/(r/r s )(l + r/r s ) 2 , 
where S c is the characteristic overdensity of the halo in terms of the critical 
density p cr i t and r s is a characteristic scale radius. 

The best fit of this function to the Coma temperature data is obtained for 
r s ~ 0.5 Mpc and S c p cr i t ~ 4 x 10~ 26 g cm~ 3 . However, this fit is formally 
unacceptable: the minimum \ 2 of 31.5 for 15 degrees of freedom can be rejected 
at greater than 99% confidence. The central temperature in the best-fit model 
is rather high ~16.5 keV but drops rapidly to ~8 keV at r s and thereafter 
continues to fall, although more gradually. The steep temperature gradient 
near the center of the cluster is inconsistent with the observed data for Coma. 
This also appears to be the case for the rich cluster Abell 2256 where the NFW 
halo function predicts a considerably steeper temperature profile in the center 
of the cluster than observed^ On the other hand, the NFW function appears 
to be an acceptable description of the optically-derived average mass profile 
of galaxy clusterso Further work on the mass distribution of Coma as well as 
other galaxy clusters is clearly needed in order to resolve this important issue. 
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